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ABSTRACT 

In this paper we have calculated the effect of Lya photons emitted by the first stars on 
the evolution of the IGM temperature. We have considered both a standard Salpeter 
IMF and a delta- function IMF for very massive stars with mass 300 Mq. We find 
that the Lya photons produced by the stellar populations considered here are able to 
heat the IGM at z ^ 25, although never above 100 K. Stars with a Salpeter IMF 
are more effective as, due to the contribution from small-mass long-living stars, they 
produce a higher Lya background. Lya heating can affect the subsequent formation 
of small mass objects by producing an entropy floor that may limit the amount of gas 
able to collapse and reduce the gas clumping. We find that the gas fraction in halos 
of mass below ~ 5 x 10^ M© is less than 50% (for the smallest masses this fraction 
drops to 1% or less) compared to a case without Lya heating. Finally, Lya photons 
heat the IGM temperature above the CMB temperature and render the 21 cm line 
from neutral hydrogen visible in emission at z ^ 15. 
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1 INTRODUCTION 

The study of the high redshift intergalactic medium (IGM) 
has recently attracted increasing attention as the new gen- 
eration of low frequency radio telescopes, e.g. L0FAl(3, Mw43 
and PAST/21cmA, promises to open an unexplored observa- 
tional window on the high-z universe. In fact, these facilities 
should be able to detect the 21 cm line associated with the 
hyperfine transition of the ground state of neutral hydro- 
gen, either in absorption or in emission against the Cos- 
mic Microwave Background (CMB) radiation (see Furlan- 
etto. Oh & Briggs 2007 for a review on the topic). In order 
for the line to be visible, the spin temperature, which reg- 
ulates the population of the levels, needs to be decoupled 
from the CMB temperature. While at z 20 collisions ef- 
ficiently decouple the two temperatures, at lower redshift 
scattering with the Lya photons emitted by the first stars 
is the most efficient process (Wouthuysen 1952; Field 1959; 
Hirata 2006; Pritchard & Furlanetto 2006). In addition, Lya 
photons could be able to heat the IGM temperature above 
the CMB temperature and render the 21 cm line visible in 
emission. The advantage of Lya photons over other heating 
sources, as e.g. x-rays (unless the spectrum is very hard), is 
that the Lya photons can travel cosmological distances and 
quickly build up a more homogeneous background. 
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The impact of Lya photon scattering on the evolution 
of the IGM temperature has been recently revised by Chen 
& Miralda-Escude (2004) and Chuzhoy & Shapiro (2007). 
Chen & Miralda-Escude (2004) included atomic thermal 
motion in their calculations, finding a heating rate several 
orders of magnitude lower than the previous estimate by 
Madau, Meiksin & Rees (1997). In addition, while "contin- 
uum" photons (with frequency between the Lya and Ly/3) 
heat the gas, "injected" photons (which cascade into the Lya 
from higher atomic resonances) cool the gas, resulting in an 
effective cooling from Lya photons at temperatures above 
10 K. Chuzhoy & Shapiro (2007) though, showed that the 
cascade which follows absorption of photons in resonances 
higher than the Lya happens via the 2s level rather than the 
2p level. Thus, the number of "injected" photons and their 
cooling efficiency are reduced compared to the estimate of 
Chen & Miralda-Escude (2004) and Lya photon scattering 
can be an efficient heating source also at temperatures higher 
than 10 K. As a consequence, Lya photons could heat the 
IGM above the CMB temperature and render the IGM vis- 
ible in 21 cm line emission. 

The IGM heating from Lya photons can also have a 
feedback effect on structure formation. In the standard cos- 
mological scenario, the first objects to collapse have typ- 
ically small masses, corresponding to virial temperature 
T^ir ^ 10* K. The formation and evolution of these struc- 
tures is heavily affected by feedback effects (for a review see 
Ciardi & Ferrara 2005). Chuzhoy & Shapiro (2007) have 
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noted that objects with virial temperatures of few thou- 
sands degrees forming in an IGM pre-heated by Lya pho- 
tons, would have a reduced post-coUapse overdensity and 
would be more vulnerable to feedback effects such as molec- 
ular hydrogen dissociation and photoevaporation. 

All this makes it worthwhile to investigate further the 
effect of Lya photons on the IGM temperature evolution. 
The paper is organized as follows. In Section 2 we esti- 
mate the Lya background expected from very high redshift 
stars; in Sections 3 and 4 we describe their effect on the 
IGM temperature evolution and feedback on early struc- 
ture formation; in Section 5 we discuss the impact on emis- 
sion/absorption of 21 cm line and in Section 6 we give our 
Conclusions. Throughout the paper we assume a ACDM cos- 
mology with = 0.76, Qm = 0.24, = 0.04, h = 0.73 
and as = 0.85. 



2 LYa PHOTON PRODUCTION 

In this Section we estimate the Lya background, Ja, ex- 
pected from the first stars. Dark matter decays and annihi- 
lations provide an additional contribution to Ja (e.g. Valdes 
et al. 2007), but this becomes negligible once the first sources 
of radiation form (Valdes, private communication). For this 
reason we limit our estimate to the stellar component. At 
each redshift, the comoving specific emissivity from such 
stars, is given by: 



(fn 



dMhdz 
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where d?n/dMhdz is the formation rate of halos of total 
mass Mh as given by the Press-Schechter formalism (Press 
& Schechter 1974), M* is the corresponding stellar mass and 
Mmin is the minimum mass of the halos. In this work we as- 
sume that stars can form in H2 cooling halos with virial 
temperature as low as 400 K and for Mmin we use the pa- 
rameterization adopted by Santos, Bromm & Kamionkowski 
(2002), i.e. Mmin{z) ~ 2.6 x lO'^ M© [(1 + z)/10]'^''^. The 
stellar mass is related to the total halo mass by: 



(2) 



where /* is the star formation efficiency and fgaa is the frac- 
tion of baryons able to collapse (see Sec. 4). In our fiducial 
case = 0.1 and fgas = 1. lv{tz^z') is the template specific 
luminosity for a stellar population of age ^/ (time elapsed 
between redshift z' and 2). We adopt the stellar spectra for 
metal-free stars computed by Schaerer (2002) and two dif- 
ferent Initial Mass Functions (IMFs): a Salpeter IMF with 
lower/upper cut-off mass of 1/100 M0, and a delta-function 
IMF for Very Massive Stars (VMS) with mass 300 Mq. For 
completeness, we have also considered a case with Pop II 
stars and a Salpeter IMF. Although the ionizing photon pro- 
duction is much reduced (by a factor of ^ 4) compared to 
metal-free stars with the same IMF, the Lya photon emis- 
sion is similar. For this reason we will limit our discussion 
to metal-free stars. 

The background intensity J^^ seen at a frequency vo by 
an observer at redshift zq is then given by: 

■^'^0(20) = r e4z)e-^^ff^''°''°'^^^dz, 




(3) 
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Figure 1. Lyo photon background (left panels) and total number 
of photons between the Lya and the Ly-limit per baryon (right) 
as a function of redshift for VMS (upper panels) and Salpeter 
(lower panels) IMF. In both panels the dotted line refers to the 
case without coupling with the feedback effect, while the solid 
(dashed) line includes the coupling and a (5„ir ~ 180 (5t,ir = 50) 
(see text for details). 



where u = i'o{i+ z) / {1+ zq) , dl / dz is the proper line element, 
Tefjiyo, zo, z) is the effective optical depth at uo of the IGM 
between redshift zo and z (see Sec 2.2 of Salvaterra & Ferrara 
2003 for a full description of the IGM modeling) and e^i^z) 
is provided by equation [T] 

In Figure [T] we plot (dotted lines) the evolution of the 
background intensity at the Lya frequency, J a, (left panels) 
and the total number of photons between the Lya and the 
Ly-limit per baryon emitted by redshift 2, Na,tot, (right). 
A case with a VMS and Salpeter IMF is shown in the up- 
per and lower panels respectively. As expected, although at 
the highest redshifts a VMS IMF produces more photons, 
a Salpeter IMF gives a higher photon density at 2 ^ 26 be- 
cause of the contribution from the small-mass long-living 
stars. 



3 IGM TEMPERATURE EVOLUTION 

In this Section we calculate the evolution of the IGM tem- 
perature, TiGM, induced by Lya photon heating. Note that 
our calculations assume that the gas is in a neutral state, i. e. 
the filling factor of ionized gas due to reionization is small. 
Just for reference, assuming an escape fraction of ionizing 
photons of 10%, we find that 10 photons per baryon are 
emitted in the redshift range (depending on other parame- 
ters adopted) 12-14 (18-20) for a Salpeter (VMS) IMF and 
8-10 for Pop II stars with a Salpeter IMF. 

TiGM is calculated with RECFAST (Seager, Sasselov & 
Scott 1999; 2000) until the first sources of radiation turn on. 
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Subsequently, we follow the method illustrated by Chuzhoy 
& Shapiro (2007). Here, we briefly summarize the relevant 
equations and we refer the reader to the original paper for 
more details. As the authors show that the presence of deu- 
terium can affect the IGM temperature (its influence in- 
creasing with decreasing redshift and increasing Na,tot), we 
have included D in our calculations. When Lya resonance 
photons are the major heating source, the H and D temper- 
ature evolve as: 
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where Ld is the D energy loss rate from collisions with H, 
TiH (no) is the number density of H (D) and Ha {His) is 
the heating rate per If (D) atom. The expressions for Ld, 
Ha and H/s can be found in the original paper and a fit 
for the energy gain/loss from "continuum" / "injected" pho- 
tons has been provided by L. Chuzhoy (Chuzhoy private 
communication). The rates Na and (the number of pho- 
tons that pass through the Lya and Ly/? resonance per If 
atom per unit time, respectively) necessary to calculate Ha 
and H/s have been derived from the model described in the 
previous Section. Our fiducial calculations assume a ratio 
of "injected" to "continuum" photons Ji/Jc ~ 0.15 and a 
deuterium abundance of [D/H]=5 x 10~^. 

In Figure [5] the temperature evolution is plotted as a 
function of redshift for a VMS (upper panel) and Salpeter 
(lower) IMF. The solid line represents the IGM temperature 
calculated with RECFAST in the absence of Lya heating, Tk- 
Th (Tb) is plotted with thin dotted (dashed) lines. As ex- 
pected, the temperatures are higher with a Salpeter IMF, 
refiecting the behavior of the Lya photon production (see 
Fig- HI- As discussed in Chuzhoy & Shapiro (2007), the most 
relevant resonance for deuterium is Ly/3. They show that D 
Ly/3 photons have a negative color temperature and as a 
consequence, they are much more efficient than Lya pho- 
tons in heating the gas. In addition, differently from H Lya 
photons, D Ly/3 photons efficiency as heaters increases with 
temperature. As a result, Td is higher than Th- We also 
find that the value of Th is independent of the presence of 
D. In the following, we will refer to Th indifferently as the 
hydrogen temperature or the IGM temperature, Tigm. 

To study the effect of the parameters adopted in our 
fiducial run, we have done additional calculations varying; 
(i) the abundance of D in the range [D/II]= 5 x 10~® — 5 x 
10~*; (ii) the ratio of "injected" to "continuum" photons in 
the range Ji/ Jc = 0.1 — 0.2; ('m^ the star formation efficiency 
/* — 0.05 — 0.2. We find that although Td strongly depends 
on [D/H] {Hp increases if [D/H] decreases), the presence 
of D does not affect Th. On the other hand, Th slightly 
increases with decreasing Ji/Jc (because the cooling from 
the "injected" photons is less effective), but the changes are 
too small to be relevant for this study. A lower (higher) value 
of the star formation efficiency induces a smaller (larger) Ja 
and gas temperature. 
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Figure 2. Temperature evolution for VMS (upper panel) and 
Salpeter (lower panel) IMF. The IGM temperature in the absence 
of hya heating, T^, is plotted (solid lines) together with Th (dot- 
ted) and Tjj (dashed). Thick (thin) lines refer to the case in which 
the feedback effect from Lya photons is included (excluded) (see 
text for details). 



4 FEEDBACK FROM LYa PHOTONS 

In this Section we will discuss the feedback effect of Lya 
heating on the formation of small mass objects. 

To study the impact of Lya feedback we follow the 
method outlined by Oh & Haiman (2003) and compare the 



IGM entropy, Kigai = Tigm/'. 



2/3 



with the entropy gen- 



erated by gravitational shock heating, Kq = Tmr /n{rmr)'^^^ , 
where uigm is IGM number density, while T„ir and n{rvir) 
are the virial temperature and the number density at the 
virial radius of a halo. The authors show that as the entropy 
parameter K — Kjgm / Kq increases, the Jeans smoothing 
effects due to the IGM entropy become more significant, the 
central pressure and density of the halos decrease, together 
with the accreted gas fraction and gas clumping (see also 
Gnedin 2000). As a consequence, the amount of gas going 
into stars is likely to decrease and the halo becomes more 
vulnerable to feedback effects such as e.g. molecular hydro- 
gen dissociation and photoevaporation. It should be noted 
that K has a minimum at r„ir, because of the higher halo 
density at r < r-^ir. Throughout this paper we will always 
be conservative and refer to the value at the virial radius. 
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Figure 3. Evolution of the entropy parameter corresponding to 
Mmin for VMS (dotted lines) and Salpeter (solid) IMF. Thick 
(thin) lines correspond to a case with Syir ~ 180 (Syir = 50). 



The value of Kq depends on the halo profile and the 
entropy parameter can be written as: 



r> _ TiGM r-2/3 
-i vir 

where Svir is the overdensity at the virial radius. For our 
fiducial runs we have adopted the value corresponding to a 
spherical top-hat collapse (Bryan & Norman 1998): 

S^ir = IStt^ + 82(0^ - 1) - 39(J7^ 

where: 



1) 



(6) 



(7) 



(8) 



For the cosmology adopted here S^ir ^ 180. 

To quantify the strength of such feedback, we plot in 
Figure [S] as an example, the redshift evolution of the en- 
tropy parameter for halos with mass Admin- The thick solid 
(dotted) Unas refer to a Salpeter (VMS) IMF. According to 
Oh & Haiman (2003) calculations (see e.g. their Fig. 6) the 
fraction of gas, fgaa, within the virial radius of a halo with 
K larger than 1 is smaller than 50% and as low as 1% for 
K ~ 10. Thus, the fraction of gas collapsing in halos with 
mass equal to Mmin can be eflectively reduced to few per- 
cents. So, to account for the eflect of Lya feedback in a self- 
consistent way, we repeat the same calculations adopting for 
/gas in equation[5]the value corresponding to K{Tigm, Mh). 
This has the effect of lowering the amount of gas that goes 
into stars and thus the intensity of the Lya background. The 
evolution of Ja and Na,tot when the coupling between the 
Lya feedback and structure formation is included, is shown 
in Figure [1] as solid lines for Svir ~ 180. As expected, both 
Ja and Na,tot are reduced by the effect of feedback. As a 
consequence, also the IGM temperature is lower when the 
coupling is included (thick lines in Fig. [2} . 
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Figure 4. Mass evolution for VMS (upper panel) and Salpeter 
(lower panel) IMF. The halo mass corresponding to T^. is plotted 
(solid lines) together with (dotted) and Mmin (dashed). 

The thick (thin) lines correspond to 6^ir ~ 180 i&vir = 50) (see 
text for details). 



Given the IGM temperature, we can derive the virial 
temperature corresponding to 7^ = 1, T^^i,. Halos with 



Tvir < T^ir ki ^ill t>e heavily affected by Lya photon feed- 
back. From T^^^ one can derive the corresponding halo 
mass, j^^. In Figure |4] the mass evolution is shown as a 
function of redshift for a VMS (upper panel) and Salpeter 
(lower) IMF. The solid lines represent the halo mass cor- 
responding to Tk. f^-^ (thick dotted lines) is plotted to- 
gether with Mmin (dashed). The feedback from Lya pho- 
tons will greatly affect the formation and evolution of halos 
with mass below Af^ , in which only about 50% of the gas 
fraction will be able to collapse. The effect is stronger for a 
Salpeter IMF. It should be noted that the Lya feedback will 
only aflect the collapse of mini-halos, while the formation of 
halos collapsing via H-line cooling will proceed unaffected. 

As mentioned above, the value of Svir depends on the 
density profile of the halo. If a Navarro, Frenk & White 
(1997) profile is adopted, the density contrast at the virial 
radius for a halo with mass Mmin varies in the range ~ 
50 — 70, depending on the redshift. As a comparison, we 
have repeated the same calculations for Svir = 50 and we 
show the evolution of J a and Na in Figure [1] with dashed 
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lines. The effect on K and ^-^ can be seen as thin lines 
in Figures [3] and [l] As expected, the impact of the Lya 
feedback is reduced compared to the fiducial case, but it 
can still be appreciable. 



5 EFFECT ON 21 CM LINE 

As already discussed in the Introduction, heating of the IGM 
from Lya photons has important consequences for the ob- 
servation of 21 cm emission line from neutral hydrogen. As 
the physics behind the emission/absorption of such line has 
been discussed extensively by several authors (for a review 
see e.g. Furlanetto, Oh & Briggs 2007), here we just write 
the relevant equations. The evolution of the spin tempera- 
ture, Ts, can be written as (Chuzhoy & Shapiro 2006): 



Ts = 



TCMB + (ya,e// + yc)Tk 



1 + yc,eff + Vc 

where Tcmb is the CMB temperature; 



■iCH + Ce + Cp) 



(9) 



(10) 



is the coupling efficiency due to collisions with H atoms, 
electrons and protons, T* = 0.068 K is the temperature 
corresponding to the energy difference between the singlet 
and triplet hyperfine levels of the ground state of neutral 
hydrogen and Aio — 2.85 x 10~^^ s~^ is the spontaneous 
emission rate. For the de-excitation rates Ch, Ce and Cp 
we have adopted the fits used by Kuhlen, Madau & Mont- 
gomery (2006; see also the original papers by Smith 1966; 
Allison & Dalgarno 1969; Liszt 2001; Zygelman 2005), i.e. 
Ch ~ nun, Ce — rie'ye and Cp — 3.2npK with nu, rie, rip 
hydrogen, electron and proton number density, k = 3.1 x 
10""T^°-3"exp(-32/rfc) cm^ s~^ and log(7e/lcm^s"^) = 
-9.607 + 0.51og(rfe)exp[-(logTfe)^ Vl800]. Finally, 

-0.3(l+z)l/^T7^/'^ 



Va.eff = 



(11) 



is the effective coupling efficiency due to Lya scattering 
which takes into account the back-reaction of the reso- 
nance on the Lya spectrum and the effect of resonant 
photons other than Lya, and ya ~ PioTi, / {AioTk) , with 
Pio ^ lO^Ja s~^ radiative de-excitation rate due to Lya 
photons. 

In the upper panel of Figure [S] the spin temperature 
is shown as a function of redshift together with Tcmb and 
Tfc (upper and lower dotted-dashed lines, respectively). The 
dotted line is Ta in the absence of a Lya background, 
while the solid (dashed) lines are for a Salpeter (VMS) 
IMF with S^ir ~ 180. The value of Tigm is hardly af- 
fected by the choice of Svir, as the difference in the Lya 
background is too small to result in an appreciable differ- 
ence in temperature. For this reason, in the following we 
will discuss only the case with S^ir ~ 180. In all cases the 
coupling with the Lya feedback is included. In the absence 
of a Lya background, Ts gets coupled to Tcmb a.t z ^ 20, 
once collisions becomes ineffective in maintaining Ts ~ Tk. 
As the intensity required for a Lya background to be ef- 
fective in decoupling the spin temperature from Tcmb is 
Ja ^ 10"^^ (1 + z) ergs cm"^ s"^ Hz"^ sr"^ {e.g. Ciardi & 
Madau 2003), Ts gets coupled to Tk as early as 2 ~ 27. 
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Figure 5. Upper panel: Evolution of the spin temperature for 
a Salpeter (solid lines) and a VMS (dashed) IMF, and in the 
absence of a Lya background (dotted) . Tcmb arid Tk are plotted 
as upper and lower dotted-dashed lines, respectively. Lower panel: 
Evolution of the differential brightness temperature. Lines are the 
same as in the upper panel. 



The differential brightness temperature, 5Ti,, between 
the CMB and a patch of neutral hydrogen with optical depth 
r and spin temperature Ts at redshift z, can be written as: 



5Tb = (1 - e 



-T\ Ts — Tcmb 
' ~z ' 



(12) 



For a gas at the mean IGM density, the optical depth is 
3Aio T*nHi(2) 



r{z) 



327rA3 Ts H{z) 



(13) 



where A is the wavelength of the transition, tihi is the num- 
ber density of neutral hydrogen and H{z) is the Hubble 
constant. 

In the lower panel of Figure [S](5Ti, is plotted for the spin 
temperature calculated above, for a Salpeter (solid lines) and 
a VMS (dashed) IMF, and in the absence of a Lya back- 
ground (dotted). From Figure [S] it is clear that, although 
Lya photons are very efficient in decoupling the spin tem- 
perature from the CMB temperature, they are not able to 
heat the IGM above Tcmb before z 15. Thus, if Lya pho- 
tons are the only heating source in the early universe and 
the reionization process is ongoing at very high redshift, we 
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expect to observe the 21 cm line in emission only when the 
IGM is substantially ionized, while at earlier times we ex- 
pect to observe it in absorption. In the absence of feedback, 
the heating by Lya photons is effective in raising the IGM 
temperature above Tcmb already at z ~ 17. 



6 DISCUSSION AND SUMMARY 

In this paper we have calculated the effect of Lya pho- 
tons emitted by the first metal-free stars on the evolution of 
the IGM temperature. We have considered both a standard 
Salpeter IMF and a delta-function IMF for very massive 
stars with mass 300 M©. Our calculations are run for a neu- 
tral gas at the mean density, thus they are strictly valid as 
long as the filling factor of ionized gas is small or in regions 
of the IGM that have not been reached by ionizing radiation. 

We find that the Lya photons produced by the stellar 
populations considered here are able to heat the IGM at 
z ^ 25, although never above ~ 100 K (thin lines in Fig. [21). 
Stars with a Salpeter IMF are more effective as, due to the 
contribution from small-mass long-living stars, they produce 
a higher Lya background. 

Lya heating can affect the subsequent formation of 
small mass objects by producing an entropy floor that may 
limit the amount of gas able to collapse. As shown in Oh 
& Haiman (2003), as the entropy parameter increases, the 
Jeans smoothing effects due to the IGM heating become 
more significant, the central pressure and density of halos 
decreases as well as the accreted gas fraction and clumping, 
with the consequence of rendering these halos more vulner- 
able to e.g. molecular hydrogen dissociation and photoheat- 
ing. We find that the gas fraction in halos of mass in the 
range 10^ — 5 x 10® Mq (depending on redshift) is less than 
50% (for the smallest masses this fraction drops to 1% or 
less) compared to a case without Lya heating. It should be 
noted that H-cooIing halos (i.e. T^ir > 10^ K) are not af- 
fected by this feedback. 

To include the effect of feedback in a self-consistent way, 
we have repeated the same calculations adopting for the 
fraction of gas able to collapse in a halo, instead of unity, 
the value induced by the feedback effect. This results in an 
overall lower amount of collapsed gas and a lower value of 
the Lya background and IGM temperature (thick lines in 
Fig. [21). To have a better understanding of the effects of Lya 
heating on structure formation one should follow the collapse 
of a halo and determine the impact on the gas density profile 
and clumping, but this is beyond the scope of this study. 

Many other feedback effects affect the formation and 
evolution of small mass halos in the early universe (see Cia- 
rdi & Ferrara 2005 for a review), ranging from photoion- 
ization/evaporation, to molecular hydrogen dissociation, to 
blowout/blowaway. Given the wide variety of feedback ef- 
fects and the fact that often they are studied in specific 
cases, it is difficult to assess the relative importance of 
Lya feedback. Differently from other feedback effects though 
(with the exception of molecular hydrogen dissociation), 
Lya photons can travel cosmological distances and quickly 
build up a fairly homogeneous background. 

Minihalo suppression from Lya feedback may also have 
an important impact on the reionization history. The lower 
star formation rate due to reduced collapse gas mass results 



in a delay of the process. On the other hand, Lya feedback 
can have the opposite effect. The halo profiles are smoothed 
out and the gas clumping factor is consequently reduced, so 
that the reionization process may be speeded up by a faster 
photoevaporation of minihalos. The final effect on the reion- 
ization history requires dedicated studies (see e.g. Shapiro, 
Iliev & Raga 2004, Ciardi et al. 2006) and is beyond the 
scope of our analysis. 

Once the IGM is heated by the Lya photons above the 
CMB temperature, the 21 cm line associated with the hy- 
perfine transition of the ground state of neutral hydrogen 
can be visible in emission against the CMB. This happens 
at 2 15, while at higher redshift, in the absence of other 
heating sources as e.g. x-rays, we expect to observe the line 
in absorption. Thus, for the detection of the emission line, 
the reionization history becomes critical as the peak of the 
emission (in terms of fluctuations of the differential bright- 
ness temperature) is expected when roughly half of the vol- 
ume is ionized but the ionized regions are not completely 
overlapped (see e.g. Ciardi & Madau 2003). To investigate 
this issue in more detail the reionization process should be 
studied together with the effect of Lya photons. It should 
be noted that metal- enriched stars would induce the same 
effect on the IGM temperature as metal-free stars with a 
Salpeter IMF. In this case though, the reionization process 
is delayed due to the much lower ionizing photon budget and 
the 21 cm line can be seen in emission over a longer redshift 
range. The same comment applies when a transition from 
an initially metal-free stellar population to metal-enriched 
stars takes place (see e.g. Schneider et al. 2006). 

To summarize the main results of this work, we find 
that: 

• the Lya background from the first stars can heat the 
IGM gas up to temperatures of ~ 100 K, depending on 
redshift and parameters adopted; 

• the associate feedback lowers the gas fraction collapsed 
in halos of mass below ~ 5 X 10'' Mq to less than 50% (for 
the smallest halos the fraction drops to 1% or less) compared 
to the case without Lya heating; 

• Lya photons heat the IGM temperature above the 
CMB temperature and render the 21 cm line from neutral 
hydrogen visible in emission a,t z 15. If Lya feedback on 
small structures is neglected emission is already possible at 
z ~ 17. 
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